finish  
/clear  
!
/filn,E354
/title,Reinforcement
*afun,deg                           ! 趨Ƕȵλ
*SET,D1,500                         ! Ͳھ  
*SET,L1,350                         ! Ͳ峤  
*SET,T,6                            ! Ͳ  
*SET,T1,6                           ! ǿ  
*SET,xd1,250                        ! ǿȦ⾶  
*SET,xd,125                         ! ǿȦھ  
*SET,xl,150                         ! ӹܳ  
*SET,xt,6                           ! ӹܺ 
*SET,P,3                            ! ʩѹֵ  
*SET,gap,0.2                        ! ǿȦͲ֮亸
*SET,ct_h1,tan(35)*T1               ! 츨Ƴߴ   
*SET,ct_h2,6
*SET,ct_h3,6
*SET,ct_h4,2
*SET,R,3                            ! ǰ뾶
qy=-P*(xd-2*xt)**2/(4*xt*(xd-xt))   ! ӹغɱ߽
qz=-P*D1**2/(4*T*(D1+T))            ! Ͳغɱ߽
/prep7  
!嵥Ԫ
et,1,solid45                        ! ʵ嵥Ԫ
et,2,TARGE170                       ! Ŀ굥Ԫ170   
et,3,CONTA173                       ! ӴԪ173
keyopt,3,5,3                        ! رշ϶
r,1,0,0,0.01,0,-gap,0               ! 趨ʵ
mp,ex,1,2.1e5                       ! 벹ǿȦϵģ
mp,prxy,1,0.3                       ! 벹ǿȦϲɱ
mp,ex,2,1.9e5                       ! ӹܲϵģ
mp,prxy,2,0.3                       ! ӹܲϲɱ
mp,ex,3,1.91e5                      ! ϵģ
mp,prxy,3,0.3                       ! ϲɱ
!Ͳӹģ
cylind,D1/2,D1/2+T,,L1,,90          ! 1/4Ͳ   
wprota,0,-90                        ! ƽת-90  
cylind,,xd1/2+T,,T1+D1/2+T,,-90     ! Ͳ,ڵõǿ⾶1 
cylind,,xd/2+xt+ct_h2,,T1+D1/2+T,,-90    !Ͳ壬ڷָǿھ9  
cylind,,xd/2+xt+ct_h4,,T1+D1/2+T,,-90    !Ͳ壬ڷָ17
bopt,numb,off   
vovlap,all                          ! 
vdele,6,11,1,1                      ! ɾ壬
cylind,xd/2,xd/2+xt,,xl+D1/2+T,,-90      ! ӹ
cylind,xd/2+xt,xd1/2,,T1+D1/2+2*T,,-90   ! Ͳ壬ڷָ23 
wpstyl,defa                              ! ָƽ
cylind,D1/2+T+gap,D1/2+2*T,,L1,,90       ! Ͳ壬ڷָǿȦԤ϶
vovlap,all                               !   
vdele,4, , ,1                            ! ɾ壬23
vdele,6,8,1,1         
vdele,15,16,1,1
vdele,19,22,1,1     
vext,80,,,0,xl,0,,,,                     ! 棬ɽӹ
wprota,0,-90                             ! ƽת-90 
cylind,,xd/2+T1+ct_h1+ct_h2,0,500,,-90   ! Ͳ壬ڷָ34  
vovlap,2,18                              ! Ͳ벹ǿȦ
vdele,5,,,1                              ! ɾ壬4
wpstyl,defa                              ! ָƽ
cylind,,D1/2+T+T1+ct_h3,,L1,,90          ! Ͳ壬ڷָ35 
vovlap,1,2                               ! Ͳӹ
vdele,7,,,1                              ! ɾ壬5
cylind,,D1/2+T-ct_h4,,L1,,90             ! Ͳ壬ڷָ18 
vovlap,1,9                               ! ͲͲ
vdele,8,,,1                              ! ɾ壬8
vext,50,52,2,0,100,0,,,,                 ! 棬Ա㴴Ͳ֮϶
vovlap,1,3,8                             ! 
vdele,9,15,6,1                           ! ɾ壬õ6
! ⺸  
askin,138,116
askin,98,72
al,15,136,13,114
al,14,137,11,115
va,3,4,7,86,6,66                         ! 湹,ͬʱõ2
!ں
askin,48,27
askin,48,102
al,16,18,78,130,96  
al,17,19,74,126,83 
va,13,14,15,8,47,51,79 
askin,101,27
vsba,12,16                               ! ú3ɵ潫ǿȦָ,õʵ
askin,65,23
vsba,13,43                               ! ú2ɵ潫Ͳָ,õʵ
! ӹܵ
afillt,30,41,R,                          ! ڽӹܵײ洴
vsba,17,58                               ! õзֽӹ
vdele,13, , ,1                           ! ɾʵ
nummrg,all                               ! ʵϢѹ
numcmp,all                               ! ʵϢ
vglue,1,4,10                             ! ճ⺸Ͳ
vglue,2,3,6,7,15,17,19                   ! ճںӹ
!ѡ
vsel,s,volu,,1,3,1
vsel,a,volu,,7,9,2
vsel,a,volu,,12,20,4
vsel,a,volu,,18,21,3
cm,weld,volu                              ! ӹ
allsel
vsel,s,volu,,23,24,1
vsel,a,volu,,5
cm,pipe,volu                              ! ӹ
allsel
vsel,s,volu,,4,8,4
cm,fill,volu                              ! ǿȦ
allsel
vsel,s,volu,,10,11,1
vsel,a,volu,,13,14,1
cm,can,volu                               ! Ͳ
!
allsel
asel,s,area,,2,5,3,1
accat,all                                 ! ͲԽǶ˲漰
esize,6                                   ! 趨ԪС
mshape,0,3D                               ! 趨Ԫ״
mat,3                                     ! 3
vsel,s,volu,,weld,1                       ! ѡк
mshkey,1                                  ! ӳ񻮷
vsweep,all                                ! ֺ
allsel
esize,12                                  ! 趨ԪС
vsel,s,volu,,pipe,1                       ! ѡнӹ
mat,2                                     ! 2
vsweep,all                                ! ֽӹ
allsel
esize,12                                  ! 趨ԪС
vsel,s,volu,,fill,1                       ! ѡвǿȦ
mat,1                                     ! 1
vsweep,all                                ! ֲǿȦ
allsel
esize,20                                  ! 趨ԪС
vsel,s,volu,,can,1                        ! ѡͲ
mat,1                                     ! 1
vsweep,all                                ! ֹ          
nummrg,all                                ! ϲԪ
numcmp,all                                ! ±

!Ӵ
allsel
asel,s,area,,19,47,28                     ! ѡͲ벹ǿȦϱ
type,3                                    ! ӴԪ
real,1                                    ! ʵ1
nsla,s,1                                  ! ѡнڵ
esln,s,0                                  ! ѡڽڵеԪ
esurf,all                                 ! ʩӽӴԪ
allsel
asel,s,area,,62                ! ѡвǿȦ±
type,2                         ! Ŀ굥Ԫ
nsla,s,1                       ! ѡнڵ
esln,s,0                       ! ѡڽڵеԪ
esurf,all                      ! ʩĿ굥Ԫ
!
/solu
!ʩԼغ
allsel                         ! ȫѡ
nsel,s,loc,z,L1                ! ѡͲ߽ڵ
sf,all,pres,qz                 ! ʩѹ
asel,s,area,,1                 ! ѡӹϱ߽
nsla,s,1                       ! ѡнڵ
sf,all,pres,qy                 ! ʩѹ
nsel,s,loc,x,0                 ! ڶԳʩӶԳԼ
dsym,symm,x
nsel,s,loc,y,0                 
dsym,symm,y
nsel,s,loc,z,0
dsym,symm,z
csys,1                         ! ϵ
nsel,s,loc,x,D1/2              ! ѡͲڱнڵ
wprota,0,-90                   ! תƽ
CSWPLA,11,1,1,1,               ! ԹƽΪ׼ûϵ
csys,11                        ! ûϵ
nsel,a,loc,x,xd/2              ! ѡнӹڱڽڵ
sf,all,pres,P                  ! ʩѹ
!趨ѡ
outres,all,all                 ! н 
time,10                        ! 趨ʱ䲽
nsub,200                       ! 趨Ӳ
nropt,full                     ! N-LѡΪFull
pred,on                        ! ʱ䲽Ԥ
lnsrch,on                      ! 
autots,on                      ! Զʱ䲽
nlgeom,on                      ! 򿪴
allsel                         ! ȫѡ
solve                          ! 
finish
!
/post1
PLNSOL, S,EQV, 2,1.0           ! 쿴ڵӦͼ
/DEVICE,VECTOR,1               ! úڰ׵ֵʾ
PLNSOL, CONT,STAT, 2,1.0       ! 쿴Ӵ״̬
PLNSOL, CONT,PRES, 2,1.0       ! 쿴Ӵѹ
